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Abstract 

We study the violation of the fluctuation-dissipation theorem in the three and four 
dimensional Gaussian Ising spin glasses using on and off equilibrium simulations. We 
have characterized numerically the function X(C) that determine the violation and 
we have studied its scaling properties. Moreover we have computed the function 
x(C) which characterize the breaking of the replica symmetry directly from equilib- 
rium simulations. The two functions are numerically equal and in this way we have 
established that the conjectured connection between the violation of fluctuation dis- 
sipation theorem in the off-equilibrium dynamics and the replica symmetry breaking 
at equilibrium holds for finite dimensional spin glasses. These results point to a spin 
glass phase with spontaneously broken replica symmetry in finite dimensional spin 
glasses. 



1 Introduction 



One of the characteristics of the disordered systems at low temperatures (and also of 
real glasses) is that its approach to equilibrium is very slow, and it is difficult to study 
equilibrium properties. Obviously in the high temperature regime there is a fast approach 
to the equilibrium. 

Due to these large time scales, the out of equilibrium regime becomes very important 
since in nature the system remains in this regime long times (minutes, days or even years). 
From the theoretical point of view it is interesting to develop a theory to describe this 
regime 

In this paper we will only discuss the low temperature phase (i.e. below the phase 
transition point of the system) and center the discussion on Ising spin glasses above their 
lower critical dimension (that lies clearly below three dimensions @). 

In the disordered case and using the Mean Field approximation (i.e. infinite range 
interactions) Cugliandolo and Kurchan have derived a generalization of the fluctuation 
dissipation theorem (FDT) that involves a new function (denoted by X) that determines 
multiplicatively (see below) the off-equilibrium regime. In the equilibrium regime X = 1 
and we recover FDT. It is possible to link this X function with the static (equilibrium) 
function x(q) (or its inverse q(x)) that appears in the replica symmetry breaking solution 
of infinite dimensional spin glasses ||. 

Unfortunately a direct check of this relation between static and dynamic in realistic 
models (like finite dimensional spin glasses) still lacks. One of the goals of this paper is to 
check this static-dynamic link in finite dimensional spin glasses. 

The crucial point of the relation between the static and dynamic is that it is possible 
to compute the complete functional form of the order parameter (the order parameter is 
a number in ordered system but it is a function, q(x), in infinite dimensional spin glass) 
using off-equilibrium simulations. Violations of the FDT relations have been reported for 
fragile glasses J3|, but in this case the corresponding equilibrium computations are still 
missing. 

On the other hand, equilibrium simulations of the three dimensional spin glasses are 
very hard and difficult . It is interesting to examine different methods than can provide 
us equilibrium information without to perform (expensive) equilibrium simulations. These 
methods exist and are based on off-equilibrium simulations (see for instance |5|, ||, [7|]). 
They have been used, for example, in the four dimensional Ising spin glass to extract the 
Edward- Anderson order parameter |J. One clear advantage is that, after a fast initial 
transient, no thermalization is needed. Another advantage is that it is possible to simulate 
large lattices and so the final results have irrelevant finite size effects. 

Following this philosophy we have computed the order parameter function^ both from 
off-equilibrium numerical simulations and from equilibrium ones and we have obtained an 
impressive agreement between both approaches that confirm the link between the static 
and dynamics in finite dimensional spin glass and provide us with off-equilibrium numer- 
ical methods to compute static quantities like the probability distribution of the overlap 
(P(g) = dx/dq) and the Edward- Anderson order parameter. 

■""We have computed directly an integrated version of the order parameter P(q), from which P(q) can 
be reobtained by double derivative. 
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We have simulated the Gaussian Ising spin glass in three and four dimensions on a 
hypercubic lattice with periodic boundary conditions. The Hamiltonian of the system is 
given by 

H = - °iJijVj • (!) 

<ij> 

By < ij > we denote the sun over nearest neighbor pairs. The are Gaussian variables 
with zero mean and unit variance. 

The plan of the paper is the following. In the next section we fix the notation and de- 
velop some analytical results. In sections three and four we show the numerical simulation 
for the three and four dimensional Ising spin glasses (respectively). Finally we present the 
conclusions. 



2 Analytical Results 

Let us fix our notations. We will study the quantity A(t) that depends on the local variables 
of our original Hamiltonian We can define the associate auto-correlation function 

C(t,t') = (A{t)A{t)) , (2) 

and the response function 

8A(t) 



R(t,t') = 



(3) 

e=0 



5e(t>) 

where we have assumed that the original Hamiltonian has been perturbed by a term 



H' = H + J e(t)A(t) . (4) 

In the dynamical framework assuming time translational invariance it is possible to derive 
the fluctuation-dissipation theorem (thereafter FDT), that reads as 

RM^m-t'f-^p- ■ (5) 

As we are interested in spin models we have chosen A(t) = <7j(t). The brackets ((■ ■ ■)) in 
eq. (|2|) imply here a double average, one over the dynamical process and a second over the 
disorder. 

The fluctuation dissipation theorem holds in the equilibrium regime, but in the early 
regimes of the dynamic we expect a breakdown of its validity. Mean Field studies H 
suggest the following modification of the FDT: 

R(t,t f )^pX(t,ty(t-t')^p-. (6) 

It has also been suggested in || § that the function X(t,t') is a function of the autocor- 
relation function: X(t,t') = X(C(t,t')). We can then write the following generalization of 
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FDT, which should hold in early times of the dynamics, the off-equilibrium fluctuation- 
dissipation relation (OFDR), that reads 

R(t,t') =/3X(C(t, t'))9{t-t')^p- . (7) 

We can relate the previous formula, eq. (|?|), with observable quantities like the magnetiza- 
tion. The magnetization in the dynamics is a function of the time and a functional of the 
magnetic field (that is itself a function of the time: h(t)) and so we can denote it m[h](t). 
Using the functional Taylor expansion we can write 



m[h](t)=m[0](t)+ I dt' 6m[h]{t) 



5h(t') 

We define the response function 

5m[h](t) 



h(t') + 0(h 2 ) 

h(t)=0 



R{t,t') 



(9) 



h(t)=0 



8h(t' 

and using the fact that in an Ising spin glass m[0](t) = 0, we obtain 

/oo 
dt' R(t,t')h{t') + 0(h 2 ) . (10) 
-oo 

Using causality we can reduce the range of the integration to (— oo,t): 

m[h](t) = f dt' R(t,t')h{t') + 0{h 2 ) . (11) 

J — oo 

This is nothing but that the linear-response theorem if we neglect the terms proportional 
to h 2 . 

By applying the OFDR we obtain the dependence of the magnetization with time in a 
generic time-dependent magnetic field (with a small strength), /i(t),Q 

m[h](t) ~ (3 J^dt' X [C(t,t')]^p-h(t') . (12) 

Now, we can perform the following experiment. We let the system to evolve in absence of 
magnetic field from t — to t — t w , and then we turn on a constant magnetic field, ho: 
h{t) = h 6(t -t w )fi Finally, with our choice of the magnetic field, we can writen 

m[h](t) ~ h p f dt' X[C(t,t')]?^p- , (13) 
Jt w at 

and by performing the change of variables u = C(t,t'), equation (O) reads 



m[h](t) ~ h P f 1 du X[u] , (14) 
Jc{t,t w ) 



2 The symbol ~ means that the equation is valid in the region where linear-response holds. 

3 Franz and Rieger jl0| used a different magnetic field function in their study of the fluctuation- 
dissipation theorem: /ifr(^) — ho9(t w — t). 

4 We ignore in our notation the fact that m[h](t) depends on t w . 
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where we have used the fact that C(t, t) = 1 (we work with Ising spins). In the equilibrium 
regime (FDT holds, X = 1) we must obtain 



m[h](t)~hoP(l-C(t,t w )) , (15) 

i.e. m[h](t)T/ho is a linear function of C(t,t w ) with slope —1. 

The link with the static is the following. In the limit t, t w — *■ oo with C(t,t w ) = q, 
X(C) — > x(q), where x(q) is given by 

x(q) = f 9 dq' P(q') , (16) 
Jo 

where P(q) is the equilibrium probability distribution of the absolute value of the overlap. 
Obviously x(q) is equal to 1 for all q > q^A, and we recover FDT. 
For future convenience, we define 

S(C) = f'dq x(q) = f'dq f dq' P(q') . (17) 
Jc JC Jo 

or equivalent ly 



d 2 C 



P(C)-- 

In the limit where X^iwe can write eq. (|14] ) as 

^^(C(M.)), (19) 
for large t w . The main goal of this paper is to test this last relation (eq. (|19D). 

3 3D Results 

The scheme of our off-equilibrium simulations has been the following. In a run without 
magnetic field we compute the autocorrelation function. We perform a second run where 
from t = until t = t w the magnetic field is zero and then for t > t w we turn on an uniform 
magnetic field of strength h Q . The starting configurations were always chosen at random 
(i.e. we quench the system suddenly from T = oo to the simulation temperature T). 

We have done a first simulation with ho = 0.1 and t w = 10 5 , with a maximum time 
of 5 • 10 6 . A second simulation was done with a smaller magnetic field, in order to control 
that linear-response works: ho = 0.05 and t w = 10 4 with the same maximum time. The 
lattice size in both cases was 64, the number of samples 4 and T = 0.7 (inside the spin 
glass phase, the critical temperature is close to 1.0 ||lT| ). 

We show in figure ([!]) the numerical results, mT/h against C(t,t w ). We have plotted 
also a straight line with slope —1 in order to control where the FDT is satisfied. 

We have also plotted the function S(C), see eq. (0), obtained at equilibrium (i.e. using 
the equilibrium probability distribution of the overlaps, P{q)) by means of a simulation of 
a 16 3 lattice using parallel tempering [|12|, [13|, ^ . We have simulated, with the help of the 
APE100 supercomputer |TJ]], 900 samples of a L = 16 lattice using the parallel tempering 



5 



method simulating 23 temperatures, from T = 1.8 down to T = 0.7 with a step of 0.05. In 
order to control the thermalization we have checked that the P(q) is completely symmetric 
in q. We have used 10 6 of sweeps (done of one Metropolis sweep and one exchange of 
temperatures) to thermalize and another 10 6 sweeps (Metropolis+Exchange) to measure 
(a detail analysis of the static of the three dimensional Ising spin glass will be presented 
elsewhere 
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Finally we have plotted two points, in the left of the figure, that are obtained with the 
infinite time extrapolation of the magnetization assuming a law 



0.18(6) and rriooT/ho 



m(t) = m c 
■ 0.46(8) in the h 



A 



0.05 run, and B 



(20) 

0.21(7) and 



with B 

m^T/ho = 0.47(4) in the h = 0.1 run. The agreement between the two Tra^/ho results 
is very good. In the statistical error there are (almost) no differences between the numerical 
curves corresponding to the two runs. 
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Figure 1: m T/h versus C with L = 64 and T = 0.7 for the three dimensional Ising spin 
glass. The curve is the function S(C) obtained from the equilibrium data. The straight 
line is the FDT prediction. We have plotted the data of the two runs: t w = 10 5 , ho — 0.1 
and t„ = 10 4 , hn = 0.05. 



From this figure we can estimate the order parameter at this temperature, that is 
precisely where the numerical curve and the straight line with —1 slope begin to be different, 
i.e. where the violation of FDT starts. We can so estimate ^ea — 0.68. We can relate this 
number with the estimate of ^ea = 0.70(2) obtained in reference [T5[] using equilibrium 
simulations. It is clear that the agreement is very good. 
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Surprisingly the S(C) curve fit very well the numerical data even in the region where 
FDT does not hold, i.e. the equilibrium distribution determines where begins the violation 
of the FDT and moreover the function x(C) is very similar to X(C) even in the very off- 
equilibrium regime, in the whole range of C . For instance S(0) = 0.45 to compare with 
the off-equilibrium data Tra^jh^ = 0.47(4). 

In this case we have been able to control down to C ~ 0.28, but with an optimal 
combination of ho and t w it should be possible to reach the region of smaller C. In any 
case the infinite time extrapolation of mT/h Q gives us the final point of the S(C) and so 
it should not be difficult to re-construct (by means of educated fits) the curve S(C) in the 
region of small C. 

This analysis implies that the Ansatz X(t,t') = X(C(t,t')) is correct in finite dimen- 
sional spin glasses and that eq. (|19|) holds in the three dimensional Ising spin glass even for 
intermediate waiting times. 



4 AD Results 

In this section we study in detail the scaling properties of the function X(T,C) and its 
dependence on the waiting time. We have used the same procedure as in the three dimen- 
sional runs. 

For the static measurements we have simulated an L = 8 lattice using the parallel 
tempering method. We have simulated 1536 samples in a range of temperatures that goes 
from T = 1.35 to 1.95 with a step of 0.05 (we remark that the transition temperature is 
1.80 ||). We have performed 10 5 sweeps (Metropolis+Exchange) to thermalize and we 
have measured, using Metropolis+Exchange, during 10 5 sweeps. This takes around one 
month of the parallel computer APE100 |)14| . We have checked that thermalization was 



achieved by analyzing the symmetry of the overlap probability distribution. From these 
simulations we have extracted the function S(C) shown in figure |2|. 

For the dynamical measurements we have performed off-equilibrium simulation using 
the same procedure that we have written in the previous section. 

We take few samples (6 in the present case) of a very large system (L = 24 and L = 32) 
such that it cannot thermalize in any computer accessible time. We have measured the 
correlation (runs without magnetic field) and the response functions of the system for 
various waiting times (t w = 2 8 , 2 11 , 2 14 , 2 17 ) verifying that for increasing t w the data of 
mT/hQ versus C(t, t w ), plotted in figure ^ collapse on a single curve loosing the dependence 
on the waiting time. We have simulated almost all the runs with h = 0.1: only in one run 
of a L = 32 lattice at T = 1.0 we have put h = 0.05. 

The clear agreement between the static and dynamical data supports (again) the cor- 
rectness of the theoretical hypothesis. Nevertheless the data for the largest waiting time 
lie a little above the static curve. We justify this discrepancy noting that in a numerical 
simulation of a relatively small volume (L = 8 in our case) the delta function for q = q E \ 
in the P(q) is replaced by a quite broad peak. This effect smoothes the cusp we expect in 
S(C) at the value C = q^A lowering the numerical curve with respect to the right one. In 
the three dimensional case the data obtained from the simulation of 16 3 lattice are very 
close to asymptotic values (by comparing, for instance, with S(C) obtained in 8 3 and 12 3 
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Figure 2: m T/h versus C with L = 32 and T = 1.35 for the four dimensional Ising spin 
glass. The curve is the function S(C) obtained from the equilibrium data. The straight 
line is the FDT prediction. Here ho = 0.1. 



lattices 0). 

Once we have verified that we can obtain information on the overlap distribution func- 
tion P(q) (measuring the linear response of a large system kept in the out of equilibrium 
regime) we have performed a systematic study in the whole frozen phase. 

We want to stress that the data from the L = 24 and the L = 32 systems coincide 
within the errors, suggesting that our results are not affected by finite size bias. Anyhow, 
we present data from both the lattice sizes. 

In figure ^| we plot the integrated response against the correlation function for different 
temperatures. The straight lines (m/ho = (1 — C)/T) represent the quasi equilibrium 
regime in which the system stays while C > (?ea- Note how the data measured in the 
regime where C < q E A collapse on a single curve independently of the temperature. 

We can understand this fact calling a hypothesis that was developed in the study of the 



P(q) in the Mean Field approximation by one of the authors (G.P.) and Toulouse [16|, 17 . 

It assumes that the order parameter q(x, T) 0, in the Mean Field theory, is a function 
of the ratio x/T for q < q E ^. This imply that we can write, in this approximation, 




for q < q E A(T) , ^ 
for q > q E A(T) , 
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Figure 3: m/ho versus C with L = 32 and different temperatures for the four dimensional 
Ising spin glass. The lines are the FDT regime: (1 — C)/T. Note how the data stay on a 
single curve when they leave the straight line (the FDT regime). Here h = 0.1. 



and, integrating x(q,T)/T, we obtain the relation between m/ho and C 

f<7EA 



rn 
h 



S(C) 



J c x(q)dq + (l-q EA )/T for C < q EA (T) 
(1 - C)/T for C > q EA (T) 



(22) 



The terms in the r.h.s. of eq. (|22|) describe the two regimes present in figure |3|: the first 
gives an expression for the curve followed by the data in the off-equilibrium regime, while 
the second is the straight line (FDT regime). 



In the next paragraphs we will show that this hypothesis [FT] also implies that the 
off-equilibrium part is independent of the temperature (i.e. in the region where C < q E \, 
[m/ho](C) is independent of the temperature). Using that the magnetic susceptibility is 
one in the spin glass phase and with the help of eq.(|2T|) it is possible, with a little algebra, 
to show that 



[1-T% EA (T))]- 



. 



(23) 



T 1 ^ y //J dT 

Now is very easy to demonstrate that the curves describing the off-equilibrium regime 
(C < Qea(T) in eq. (|22D ) do not depend on the temperature. By deriving the curve expres- 
sion with respect to T we obtain 



d 


"m" 


d 


~s(cy 


dT 


.ho. 


~ dT 


T 



x(q EA (T)) 



dq EA (T) 1 dq EA (T) 



Qea(T) 



dT 



T dT 



J^2 



(24) 
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where in the last equality we have made use of eq.(^). So we have verified that the 
first expression in eq.(E3) does not depend on T. We finally write that for C — > q EA the 



hypothesis [16, O] implies 



S(C) ~ VT^C . (25) 

At this point we have seen that Mean Field predicts qualitatively the behavior plotted 
in figure |^ for a finite dimensional spin glass. Now we will examine quantitatively the data 
of figure H 

For C < q EA we have seen (figure |3|) that the numerical data can be approximated by 
a power law of the variable 1 — C 



h 



j TA(1 - C) B for C < q EA (T) , 
~ 1 1 - C for C> q EA (T) , 



(26) 



with A ~ 0.52 and B ~ 0.41 (not very far from the Mean Field behavior, (1 — C) 1 / 2 ). 
Multiplying both sides of the previous expression by T -1 ^ 1- ^ we have 



mT^^. T~ — A(1-C) B = A 

-T l-B 



(l-C)T^j for C<q EA (T) , 
h o ' { T~^(l-C) = (l-C)T-* forC>g E A(T), 



(27) 



where we have introduced (j) = 1/(1 — B) ~ 1.7 for convenience. Doing so we can rescale 
the data for all the temperatures on a single curve like the one shown in figure [|. 

The good scaling of data (figure |j) obtained with different magnetic fields is a con- 
firmation that we are working in the linear response regime. It should be noted also the 
absence of different finite size effects for the lattices we have considered (24 4 and 32 4 ). 



5 Conclusions 

In this paper we have found that the violation of the fluctuation-dissipation theorem in 
finite dimensional spin glasses follows the lines of the violation of the theorem in Mean 
Field models. 

We have also found that the function that determine the violation is given, since inter- 
mediate waiting times, by the double integral of the probability distribution of the overlap 
calculated at equilibrium. 

This fact gives us a further confirmation that the Ansatz used in references || are 
correct even in finite dimensional models (i.e. X depends only on C, as it was established 
in reference |10j)- We have also obtained that the violation of the theorem is given by the 
static (i.e. we can express X{C) as a function of static quantities). 

Moreover we have seen that by controlling the scaling of the waiting times it is possible 
to construct the X(C) curve without doing equilibrium simulations. Also these curves 
provide us an useful and precise method to compute the Edward- Anderson order parameter. 

The form of the X(C) function is very different from that of droplet (a distinguish 
ferromagnet), ferromagnetic and one step replica symmetry breaking systems ||, and so 
we have obtained another evidence that the finite dimensional Ising spin glasses cannot be 
described by the droplet model. 
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Figure 4: (mT/ho) versus (1 — C) T'^ with <fi = 1.7. Note that in the plot we have 
included data measured on different lattices and in presence of different magnetic fields. In 
the FDT regime (left part of the figure) the factor has no effect because in this region 
mT/ho = 1 — C. The off-equilibrium regime (right part of the figure) follows a power law 
with power B = 0.41. 
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Finally we have studied the scaling properties of X(C) finding that it is possible pa- 
rameterize it using static Mean Field analytical results. It gives us a further evidence of 
spontaneously broken replica symmetry (infinite steps of replica symmetry breaking). 
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